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Gel electrophoresis is a powerful experimental method to probe the 
topology of DNA and other biopolymers. While there is a large 
body of experimental work which allows us to accurately separate 
different topoisomers of a molecule, a full theoretical understand¬ 
ing of these experiments has not yet been achieved. Here we show 
that the mobility of DNA knots depends crucially and subtly on 
the physical properties of the gel, and in particular on the pres¬ 
ence of dangling ends. The topological interactions between these 
and DNA molecules can be described in terms of an “entanglement 
number”, and yield a non-monotonic mobility at moderate fields. 
Consequently, in two-dimensional electrophoresis, gel bands display 
a characteristic arc pattern; this turns into a straight line when the 
density of dangling ends vanishes. We also provide a novel frame¬ 
work to accurately predict the shape of such arcs as a function of 
molecule length and topological complexity, which may be used to 
inform future experiments. 

DNA Knots | Gel Electrophoresis | Topology 

T opology plays a key role in the biophysics of DNA, and is 
intimately related to its functioning. For instance, tran¬ 
scription of a gene redistributes twist locally to create what is 
known as supercoiling, while catenanes or knots can prevent 
cell division, hence they need to be quickly and accurately re¬ 
moved by specialised enzymes known as topoisomerases. But 
how can one establish experimentally the topological state of 
a given DNA molecule? By far the most successful and widely 
used technique to do so is gel electrophoresis [1, 2]. This 
method exploits the empirical observation that the mobility of 
a charged DNA molecule under an electric field depends on its 
size, shape and topology [2]. Gel electrophoresis is so reliable 
that it can be used, for instance, to map replication origins and 
stalled replication forks [3], to separate plasmids with different 
amount of supercoiling [3, 4], and to identify DNA knots [5, 6]. 
The most widely employed variant of this technique nowadays 
is two-dimensional gel electrophoresis, where a DNA molecule 
is subjected to a sequence of two fields, applied along orthogo¬ 
nal directions [2]. The two runs are characterised by different 
field strengths, and sometimes also gel concentrations [4]; with 
suitable choices, the joint responses leads to increased sensi¬ 
tivity. 

While gel electrophoresis is used very often, and is extremely 
well characterised empirically, there is still no comprehensive 
theory to quantitatively understand, or predict, what results 
will be observed in a particular experiment. Some aspects are 
reasonably well established. For instance, it is now widely ac¬ 
cepted that the physics of the size-dependent migration of lin¬ 
ear polymers can be explained by the theory of biased polymer 
reptation [7, 8, 9, 10, 11, 12]. Likewise, the behaviour of, for 
example, nicked, torsionally relaxed, DNA knots in a sparse 
gel and under a weak field is analogous to that of molecules 
sedimenting under gravity [13, 14]. The terminal velocity can 
be estimated via a balance between the applied force and the 
frictional opposing force, which is proportional to the average 
size of the molecule: as a result more complex knots, which are 
smaller, move faster under the field. However, the mechanisms 
regulating the electrophoretic mobility of DNA knots at inter¬ 


mediate fields, and in more concentrated agar gels, is much 
less understood [4, 12, 15]. Here, experiments suggest that 
the mobility of knots is usually a non-monotonic function of 
the knot complexity, or, more precisely, of their average cross¬ 
ing number [16, 5] (ACN): initially knots move more slowly as 
their ACN increases, while, past a critical ACN, more complex 
knots move faster. The combination of the responses to ex¬ 
ternal fields directed along two perpendicular directions leads 
to a characteristic electrophoretic arc which allows to separate 
the first simple knots more clearly in a 2D slab [17, 18, 6, 4]. 
To our knowledge, there is currently no theoretical framework 
that quantitatively explains the non-monotonic behaviour at 
intermediate or large fields and the consequent formation of 
arc patterns. 

To address this issue, here we present large scale Brown¬ 
ian dynamics simulations of knotted DNA chains migrating 
through a gel, and subjected to a sequence of fields of dif¬ 
ferent strength and direction, as in two-dimensional gel elec¬ 
trophoresis experiments (see cartoon in Fig. 1(a)). We model 
the gel as an imperfect cubic mesh [19], where some of the 
bonds have been cut (see Methods) to simulate the presence 
of open strands, or dangling ends, which have been observed in 
physical agarose gels [20, 21, 22, 23, 24]. Our results confirm 
the linear relation of the electrophoretic mobility with average 
crossing number for the first simple knots (we study ACN up 
to 12) in a sparse gel and under a weak field. However, our 
simulations also suggest that, due to a non-negligible proba¬ 
bility of forming “impalements” where a dangling end of the 


Significance 

Gel electrophoresis is a ubiquitous biophysical technique. It con¬ 
sists in dragging charged biopolymers through a porous gel, by 
applying an electric field. Because the migration speed de¬ 
pends on topology, this method can be used to classify DNA 
knots. Currently, electrophoresis largely relies on empirical ob¬ 
servations, while its theoretical understanding is limited. For 
instance, no theory can explain why knot mobility under strong 
fields depends non-monotonically on complexity. Our study re¬ 
veals a possible reason: complex knots have a smaller size, hence 
move more quickly through the gel pores; however, once they 
get trapped by the dangling ends present in the gel, the resulting 
entanglements are slower to resolve. Our results can improve 
the design of future electrophoresis experiments. 
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gel pierces a knot, the response of the chain to stronger fields 
is different. We suggest that, in this regime, the sole radius 
of gyration is not enough to explain the observed dynamics, 
and we introduce an average “entanglement number” which 
increases with the ACN and provides a measure of the likeli¬ 
hood of forming an impalement. The time needed by a knot 
to disentangle from the gel increases with its average entan¬ 
glement number, or knot complexity, and this slows down the 
motion, thus competing with the Stokes friction which leads to 
an increase of mobility with ACN. As a result of this contest, 
one typically gets a non-monotonic behaviour of the terminal 
speed with ACN, and an electrophoretic arc in 2 dimensions. 

We also provide a simple model, based on a mapping be¬ 
tween the DNA knot dynamics and a biased continuous time 
random walks, which faithfully reproduces our Brownian dy¬ 
namics simulations starting from a minimal number of as¬ 
sumption. This approach can then be used to predict how 
the shape of the electrophoretic arc should depend on system 
parameters such as the average lattice spacing (pore size) of 
the gel and the contour length of DNA knots; as we shall 
see, the predicted trends are in agreement with existing two- 
dimensional electrophoresis data. This constitutes, to our 
knowledge, the first example of a quantitative prediction of 
two-dimensional electrophoretic diagrams, hence we suggest 
that the approach we present could potentially lead to even 
more accurate and targeted experiments to separate topoiso- 
mers in DNA or other polymers. 



Fig. 1. (a) Snapshot (to scale) of the model gel with some examples of knotted 

configurations. Note that to model a physical gel, a simple cubic structure is ran¬ 
domly cut to create dangling ends. (b) Equilibri urn configurations of some of the 
knots considered; it can be readily seen that the size tends to be smaller as the knot 
becomes more complex. The knots pictured in (a) and (b) are: trefoil (3i) in red, 
figure of eight (4i) in orange, pentafoil (5i) in yellow. Stevedore’s (6i) in green and 
“nonafoil” (9i) in grey. 


DNA knots form an electrophoretic arc only in irregular gels. 

The system we studied, sketched in Fig. 1 (see also Materials 
and Methods and SI), consists of 10 nicked, i.e. torsionally 
relaxed, DNA loops of ~ 3.7 kilo-base pairs (kbp) within a 
model agarose gel with dangling ends. These loops are either 
unknotted, or form one of the first few simple knots (with up 
to 9 crossing in their minimal projection [25]). The loops were 
first equilibrated within the gel (see SI), and then subjected 
to an in silico gel electrophoresis process where a weak elec¬ 
tric field is first applied (< 50 Vjcm) along the vertical (z) 
direction, followed by a stronger field (> 150 Vj cm) along a 
transversal, say y, direction. We refer to these two fields as 
“weak” and “stronger”, or “moderate”, in what follows. The 
complete equations of motions and force fields used in our 
Brownian dynamics simulations are detailed in the SI. 

By monitoring the trajectories of the knots through the gel 
we computed the average speed of their center of mass along 
each of the field directions (see SI, and Fig.2a-b). As expected 
the mobility along the direction of the weak field increases 
with the topological complexity of the configurations. Along 
the direction of the moderate field, however, the mobility of 
the knots displays a non-monotonic behaviour. In particular, 
the unknot now moves faster than either the trefoil or the 4i 
twist knot, and has an average speed similar to the 5i knot. 
This non monotonic behaviour of the knot mobility, as a func¬ 
tion of the ACN, was previously observed in typical experi¬ 
ments with torsionally relaxed DNA knots [17, 18, 6, 26, 4]. 
It is worth noticing that, within the 5-crossings family the 5i 
torus knot moves more slowly than the 52 twist knots. This 
is similar to what observed in the weak field case (although 
much less enhanced) but different from what is observed in 
experiments of sedimentation [14]. 

To better compare our findings with experiments we report, 
in Fig. 2c the spatial distribution of the knots as Gaussians 
centred in (vz{lC)tz-,Vy(lC)ty) where Vz(lC) and Vy(lC) are the 
velocities along z (weak field direction) and y (stronger field 
direction) of knot K, while ty and tz the electrophoretic run¬ 
times. The width of the Gaussians is set to be proportional to 
the standard deviation of the velocities. The resulting spots 
can be seen as the in silico analogue of the ones observed in 
gel electrophoretic experiments. Note that the combination of 
a monotonic behaviour along the weak field direction with a 
non-monotonic one along the stronger field, gives rise to the 
arc shape distribution of the spots characteristic of 2D elec¬ 
trophoresis experiments run either on knotted configurations 
or on or supercoiled plasmids [4] (see Fig. 2d). 

It is interesting to ask whether one can observe the elec¬ 
trophoretic arc also in simulations where the gel is a regular 
cubic mesh, i.e. a mesh with no dangling ends, as this has 
been so far the typical way to model an agarose gel [15]. Re¬ 
markably, unlike the case of gel with dangling ends, also called 
“irregular” hereafter, no example of non-monotonic behaviour 
of the knot mobility is found for regular gels (for comparison 
see Fig. 2c) and Fig. SI in SI). This result is in line with 
previous simulations based on lattice knots in regular gels [15] 
and persists for different field strengths (1.25 — 600 Vj cm) and 
gel pore sizes (200 — 500 nm) (see SI). ^ Hence our simulations 
strongly suggests that the causes for the non-monotonic be¬ 
haviour, observed in the case of irregular gels, are to be sought 
in the interaction between the knots and the gel dangling ends. 

This conjecture is also supported by the fact that lin¬ 
ear (open) DNA samples, are frequently observed to migrate 


^We note that physical gels also have an inhomogeneous pore size; while considering this aspect will 
affect our results quantitatively, the common understanding is that the knot speed should depend 
monotonically on size [1]. This is qualitatively different from the case of dangling ends, where the 
gel-polymer interactions strongly depend on topological complexity as well. 
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(b) acn 




Fig. 2. In silico 2D gel electrophoresis for knotted DNA loops, (a) Average velocity of different types of knotted polymers along the direction of the weak electric 
field: Ei = fijqb — 50 Vjcm. The dotted line indicates the linear increase with the ACN. (b) Average velocity of different types of knotted polymers along the 
direction of the moderate field E 2 = f^/qb — 150 Vjcm. (c) 2D reconstruction of the spatial distribution of the knots as Gaussian probability distributions centred 
in {Vz{E)tz,Vy{)C)ty) where Vz(JC) and Vy{lC) are the average speeds along, respectively, the weak and stronger field direction of knot type K, as found by Brownian 
dynamics simulations. The electrophoretic run-times correspond to tz = 0.8 10^^ TBr — 5 min and ty = 2.6 10^® tbt — 16 min. The spread of each spot has been 
estimated by looking at the standard deviation of Vz{lC) and Vy{lC). (Examples of single trajectories along the 2 : and y directions are reported in Fig. SI), (d) Outcome 
of a 2D gel electrophoresis experiment performed on P4 viral DNA (10 kbp) at 0.4% agarose concentration (reproduced from Ref. [6] with permission). 


faster than covalently closed (unknotted) ones, in gel elec¬ 
trophoresis experiments performed in both strong and weak 
fields [17, 6, 26, 4]. This is in line with the outcome of a re¬ 
cent computer experiment probing the dynamics of linear and 
unknotted circular molecules through an irregular gel [19]. 


DNA samples become severely entangled with the dangling 
ends. Having established that the presence of dangling ends 
in gels affect severely the transport properties of the knotted 
DNA loops under moderate electric fields, it is natural to look 
at the possible mechanisms ruling this phenomenon. 

The typical trajectories and average extension of some knot¬ 
ted loops, as they move through a regular model gel and a gel 
with dangling ends are markedly different at moderate fields 
(see S.I. Fig. SI). In a regular gel, knots respond to the field, 
by shrinking their size so as to channel through the pores of 
the gel more efficiently. This mechanism, also known as “chan¬ 
nelling” , for which polymers squeeze through the gel pores, has 
already been observed in previous works [27, 28], and it was 
previously conjectured to play a role in the non-monotonic 
separation of DNA knots in gels, as more complex knots could 
have a different ability in deforming their overall shape when 
squeezing through the pores [4]. On the other hand, as dis¬ 
cussed in the preceding section, we find that this behaviour is 
not sufficient to explain the electrophoretic arc, as for regular 
gels we always observe a monotonic separation of the knots as 
a function of the ACN (see S.I., Fig. SI). On the contrary, in 
the case of irregular gels, knotted loops are much more prone 
to entangle with one (or more) dangling ends (see insets of 
Fig. 2c for some examples). These entangled states (or “im¬ 
palements” ) require some time to be unraveled and this is the 
reason of the anomalously long pauses observed in the knot 
trajectories (see in particular SI Fig. SI). Clearly, as the 
DNA gets longer, “impalements”, which can either be par¬ 


allel or perpendicular with respect the direction of the field, 
become progressively more likely. As a matter of fact this 
could be one of the reasons why it is in practice unfeasible to 
perform efficient gel electrophoresis experiments with circular 
DNA longer than 10 kbp [29]: at these sizes impalements are 
so frequent that they may cause DNA breakage. 

In analogy with the phenomenon of threading, which slows 
down the dynamics of unknotted loops either in a melt or in 
a gel [30, 31, 19], and that of “crawling” of knots around ob¬ 
stacles [32], it is reasonable to expect that more complicated 
knots will take longer to disentangle themselves from an im¬ 
palement. We argue that this mechanism, when competing 
with the reduced Stokes drag of more complex knots in gels, 
is ultimately responsible for producing a non-monotonic de¬ 
pendence as a function of their complexity, i.e. their ACN. 

More complex knots have smaller size but larger ‘‘entangle¬ 
ment number”. Given that “impalement” events are key fac¬ 
tors in determining the mobility of DNA knots within gels 
with dangling ends, it is important to find a way to define and 
measure this entanglement. Impalement may occur with dan¬ 
gling ends oriented along several directions (see examples in 
Fig. 2(c)) but it is reasonable to expect that all these events 
involve a similar mechanism in which, i.e. one dangling end 
“pierces through” the knot. 

To quantify the degree of knot-gel entanglement, we con¬ 
sider an equilibrated knot configuration in the gel, and project 
it on the plane perpendicular to the field direction. We then 
choose randomly a base point P, at a distance from the pro¬ 
jection plane that is much bigger than the radius of gyration 
of the projected configuration, Starting from P we 

draw an arc which pierces only once the projection plane at a 
point, Q, chosen randomly, with uniform probability within a 
disk of radius and centred in the center of mass of the 

projected configuration. The arc and the plane define a semis- 
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Fig. 3. (a) Sketch of the procedure used to define the piercing, or entanglement number: for a given projection of the configuration, the crossings define a set of regions 

pvoj 

whose intersections with the disk of radius Rg are highlighted in blue. Starting from different points P far away from the projections closed paths (2 in the example) that 
pierce once the disk at different locations Q are built at random. The absolute value of the linking number is then computed between the knot configuration and each closed 
path. The average over the set of closed paths is finally taken: this is the average entanglement number, (b) and (c) show the results of the procedure described in panel 
(a) for a configuration with knot type 9i (c) and for an unknotted configuration (b). The regions are coloured according to the computed absolute value of the linking 
number (see colour map at the left). Note that for the 9i case there are regions of high (3) \Lk\, which are more prone to become entangled with the dangling ends of the 
gel. (d) shows the average entanglement number (tt) and the mean squared radius of gyration divided by (//2)^ {I is the gel pore size) for different knot types classified in 
terms of ACN. 


pace and we close the path with a second arc connecting Q 
and P and living in the other semispace (see Fig. 3a). To as¬ 
sess whether this circular path interacts topologically with the 
knotted configuration, we compute the absolute value of the 
linking number, \Lk\, between the circular path and the knot. 
Fig. 3b and Fig. 3c show the result of this procedure when 
applied to two different knotted loops. By averaging \Lk\ over 
several circular paths with different Q, and over different knot 
conhgurations, we dehne the “average entanglement number”, 
(tt) (AEN) as the measure of the degree of entanglement be¬ 
tween the knotted loops and the surrounding gel. 

From Fig. 3d we see that (tt) grows approximately linearly 
with ACN: as one would expect, more complex knots can, on 
average, become more entangled with the surrounding irregu¬ 
lar gel. It is interesting to compare this behaviour with that of 
the mean squared radius of gyration normalised with respect 
to half the gel pore size (see Fig. 3d): unlike the average en¬ 
tanglement number, the average extension of the loop is, to a 
good approximation, inversely proportional to the ACN, i.e. 
to the knot complexity. This corresponds to the well-known 
fact that, for a given loop contour length, more complex knots 
are on average less extended [5, 14] (see also the equilibrium 
conhgurations in Fig. 1). 

The plots in Fig. 3d suggest a possible interpretation of the 
non monotonic mobility of the knots in irregular gels based 
on the interplay between the average size and the degree of 
entanglement with the gel. On one hand more complex knots, 
being smaller in size, experience less frequent collision with 
the gel and hence should travel more easily through it: this is 
just another variant of the Stokes friction argument discussed 


previously. On the other hand, once knot-gel collisions occur, 
more complex knots experience a more intricate entanglement 
with the gel (higher values of AEN are more probable) that 
will take longer to be unravelled [5, 13, 15, 32]. 

The above argument suggests the existence of two time 
scales in the process: one is the time r/ between two suc¬ 
cessive knot-gel collisions yielding a local entanglement; the 
other, Tdis, is the time needed by the knotted loop to fully 
disentangle from the impalement. The time scale r/ increases 
as the knot average size decreases and hence increases with 
knot complexity (ACN). In other words more complex knots 
experience, on average, less collisions with the gel than their 
simpler counterpart. The second timescale, Tdis, is instead an 
increasing function of (tt) (see Eig. 3(d)) and hence of the 
knot complexity (measured in terms of ACN). According to 
this picture, the slowest topoisomer in an irregular gel with a 
given lattice spacing will be the one with the “best” compro¬ 
mise between a high rate of collisions, and a sufficiently high 
value (tt). 

To investigate more quantitatively the dependence of r/ 
and Tdis on the knot type (ACN), we analyse the trajecto¬ 
ries of the knotted loops in the gel by computing: (i) the 
average number of times a knot arrests its motion in the gel 
(entanglement event), (rie), and (ii) the distribution of the 
duration of these entanglement events. As specified in the 
SI, the duration of the entanglement events can be identified 
as the time intervals where the spatial position of the centre 
of mass of the conhguration deviates significantly from the 
expected collision-free field-driven linear motion with speed 
Vfree — Fz,y/^C, — fz,y/C- 
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Fig. 4. (a) Average number of events in which the knot is entangled with the 

surrounding gel (entanglement events) as a function of ACN. (b) Average disentan¬ 
glement time as a function of ACN. In these estimates only entanglement events with 
duration shorter than 200 Tsr ai'e considered. 

As reported in the SI (Fig. S 2 ), the average fraction of 
time in which the knot is trapped, Twlrtot {rtot is the time of 
the full trajectory), is a non-monotonic function of the ACN, 
in line with the result on the mobility under moderate field 
(Fig. 2b). In Fig. 3(a) we show that the average number of 
entanglement events {rie) decreases with the knot complexity 
(i.e. ACN). On the other hand the distribution of the du¬ 
ration of these events displays an intriguing bimodal shape, 
with two peaks occurring respectively at short and very long 
times (see SI Fig. S3). The peak at long times can be in¬ 
terpreted as the signature of head-on collisions with the gel, 
where the dangling ends involved are opposite to the direc¬ 
tion of the knot motion. The resulting entanglement is, in 
this case, very difficult to unravel especially in presence of a 
strong electric field (see for instance the inset of Fig. 2 for 
the unknot). Nonetheless, either if we exclude or not from the 
statistic the entanglement events corresponding to the peak at 
long times the characteristic disentanglement time Tdis turns 
out to increase (linearly) with the knot complexity, i.e. with 
ACN (see Fig. 4(b)). 

We assume that this bi-modal shape is due to a shift in 
the energy barrier that the knots have to overcome in order 
to disentangle from the dangling ends. In particular, one can 
think to this process as an Arrhenius process, where the en¬ 
ergy barrier is a function of the length (projected along the 
field direction) of the dangling end, the knot complexity and, 
more importantly, the magnitude of the external field. When 
this is too strong, disentanglement events are very rare, and 
all knots will end up being permanently entangled with the gel 
structure; on the other hand, when this is too weak, the typ¬ 
ical disentanglement time is very short, and the dependence 
of Tf as a function of the ACN dominates the motion of the 
polymers, re-establishing the usual linear relationship. 


A random walk model with topology-dependent rates cap¬ 
tures the observed non-monotonic behaviour. As shown in 
previous sections, 2D electrophoresis experiments and Brown¬ 
ian dynamics simulations of knotted loops in irregular gels are 
in qualitative agreement under many aspects. In this section 
we propose a simple model that reproduces the main findings 
of the simulations and furnishes a simple but accurate way 
to predict the arc shapes of the experimental patterns as a 
function both of the knot complexity and of the loop contour 
length. 

In this model, we describe the knotted loop moving within 
the irregular gel as a biased random walk on a ID lattice i.e. a 
random walk that moves to the right (direction of the external 
field) unless it is trapped into an entangled state (due to im¬ 
palement) with probability Ae(/C) = (see SI for more 

details). Once in the entangled state, the walker has to wait a 
given amount of time that is picked randomly from a bimodal 
distribution consisting of an exponential decay, modelling the 
short time disentanglement, and a smaller probability peak at 
large times, describing the long disentanglement time from a 
head-on collision (see SI for the details). In this simple descrip¬ 
tion the only relevant parameters are the hitting rate and the 
parameters characterising the bimodal distribution of waiting 
(i.e. disentanglement) times. 

Once the values of these parameters are set to reproduce the 
data reported in Fig. 3 (see also SI), the model can be used 
to predict the mobility of the electrophoretic arc as a func¬ 
tion of ACN. As shown in Fig. 5 this procedure reproduces 
with remarkably good agreement the simulation data and, in 
particular, it captures the physical mechanism leading to the 
non-monotonic mobility at moderate field. Note that, as the 
random walker solely moves to the right, the field strength 
enters into the model only through the waiting times and the 
hitting rates. 

More importantly, once the parameter values of the biased 
random walk model are set for a given pore size of the gel, 
/i, their values for a different pore size, I 2 , can be estimated 
from general arguments (see SI). We can therefore use this 
simplified model to predict the moderate field mobility and 
the shape of the electrophoretic arcs of DNA knots in gels of 
variable pore size, e.g. tuned via agarose concentration [33] or 
nano-wire growth cycle [34]. 

The plots presented in Fig. 5(b) suggest that tighter gels 
give rise to more curved (or deeper) arcs where the slowest 
knot has a higher ACN with respect to sparser gels. More¬ 
over, since the entanglement rate Ae and the disentanglement 
time Tdis (both of these relative to the same quantities for the 
unknot) should depend only on the ratio between the knot 
extension and the gel pore size /, a similar trend should be 
observed also by increasing the DNA loop contour length by 
keeping fixed I (see SI). This is in qualitative agreement with 
experiments, as electrophoretic arcs are straighter for shorter 
DNA molecules (Fig. 5(d)). A further quantitative prediction 
we can draw from our arguments is that the relative position 
of the three 6 -crossing knots can be controlled by tuning the 
pore size (Fig. 5) of the gel. Indeed the size of the pores deter¬ 
mines whether the 61 Stevedore’s knot is to the left or to the 
right of the minimum of the mobility curve: in the former case 
61 will moves faster in the gel than the 62 and 63 knots (which, 
having higher ACN have also higher AEN), while in the latter 
case it will move more slowly. This detailed prediction could 
be tested in future electrophoresis experiments with knotted 
DNA loops moving within different gels. 

Conclusions. We have studied the role of topology in the gel 
electrophoretic mobility of DNA knots by means of Brownian 
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Fig. 5. (a) Average speed along the direction of the moderate field from Fig. 2(b). The dashed line is obtained from the biased continuous random walk model, and 

corresponds to the (shifted and rescaled) red curve in (b). (b) Average relative separation (in units of lattice spacing over time) of the knots as a function of the ACN for 
different parameters, as predicted by the continuous random walk model. The grey dashed line in (a) is obtained by shifting the red curve in (b) by the value of (t’2/)(0i) 
and rescaling it by the free velocity Vf^ee- (c) Reconstruction of a 2D gel electrophoresis experiment from the data in (b) and zoom over the relative position of the family 
of 6-crossings knots for two cases in which the minimum of the arc is at their left and their right, (d) Outcomes of a 2D gel electrophoresis experiments performed on P4 
viral DNA with different lengths, respectively 4.7 (black) and 10 (white) kbp at equal agarose concentration (0.4 %) (reproduced from Ref. [26] with permission). 


dynamics simulations and a minimal model of biased random 
walk. We showed that, when the knots are driven through a 
physical gel i.e. possessing dangling ends, the knots’ mobility, 
as a function of their average crossing number (ACN), depends 
on the strength of the external field. 

At weak fields we recover the well-known linear relationship 
between migrating speed and knot-type; at stronger fields, we 
observe instead a non-monotonic behaviour. We argue that 
this puzzling feature, routinely observed in experiments but 
yet not fully explained, can be better understood by taking 
into account the topological interactions, or entanglements, of 
the knots with the irregularities of the surrounding gel. While 
more complex knots assume more compact configurations, and 
hence smaller Stokes friction than simpler knots, they also ex¬ 
perience more complex entanglements with the gel and hence 
longer disentanglement times. These two competing effects 
give rise to the non-monotonic speed of the knots observed 
in the experiments, a feature that, remarkably, is absent for 
knotted loops moving in a regular gel {i.e. no dangling ends). 
While most of our simulations were performed with a rigid gel, 
we tested that the results are qualitatively unchanged for gels 
with flexible dangling ends (see Materials and Methods, and 
SI). 

We also propose a model which describe the motion of knot¬ 
ted DNA loops as a biased continuous time random walks. 


This model, although minimal, by focusing on the competition 
between Stokes friction and topological entanglements high¬ 
lighted by the simulations, is able to predict the shape of elec¬ 
trophoretic patterns of DNA knots of different contour length 
observed in gels with tunable physical properties. In particu¬ 
lar we predict that, by changing the ratio between the radius 
of gyration of the unknot and the gel pore size, 2D gel elec¬ 
trophoresis experiments should lead to deep electrophoretic 
arcs for tight gels (or long knots), and shallow ones for sparse 
gels (or short knots). 

We hope that our results will prompt further experimen¬ 
tal and numerical verification on the role of topology in the 
anomalous electrophoretic mobility of knotted polymers and, 
consequently, suggest new and more accurate setups to sepa¬ 
rate biopolymers of different topology. Lastly, it is likely that 
a similar competition between loop size and loop-gel interac¬ 
tions can be responsible of the observed electrophoretic arcs 
in DNA molecules with different density of supercoiling; we 
plan to explore this related situation, as well as the case of 
composite knots, in the future. 

Materials and Methods 

Double-stranded (ds) and nicked, i.e. torsionally relaxed, DNA knots are modelled 
as closed and knotted semi-flexible bead-spring chains [35], with beads of diameter 
a = 2.5 nm, which reflects the thickness of hydrated B-DNA near physiological 
conditions [36]. The persistence length is set to Ip = 20a = 50 nm, and the 
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chosen contour length Lc = 512(7 corresponds to DNA loops of length ~ 4 kbp 
~ 1.3 fi m. 

The gel is modelled as an imperfect and rigid cubic mesh, with lattice spacing 
I = 80(7 ~ 200 nm compatible with the average pore size of agarose gels at 5% 
and artificial gels made of solid nano-wires [33, 34] (for more details on the model 
and comparison with the case of flexible dangling ends, see SI). The irregularities, or 
dangling ends, of the gel are created by starting from a regular cubic mesh and then 
by halving some of the edges randomly, with probability p = 0.4. Although this 
probability is chosen arbitrarily, it is possible to map it to a real value of “disorder” 
found in an agarose gel at a given concentration by comparing the mobility of linear 
and ring polymers running through it, similarly to what was done in Ref. [19]. The 
edges of the mesh are discretised with beads of size ag = 10a ~ 25 nm which 
is compatible with the observed diameter of agarose bundles [33, 37]. For simplicity 
and computational efficiency, we model the gel as a static mesh, meaning that the 
mesh structure is not deforming under either thermal of mechanical strains. This is 
an approximation for an agarose gel, whose bundles are generally found, at the con¬ 
centrations used in gel electrophoresis, to be made of tens of fibers whose persistence 
length has been observed to be around 2 — 10 nm [37]. In light of this, a conservative 
estimation of the persistence length of an agarose bundle is comparable with that of 
DNA, i.e. Ip ~ 50 nm (this assumes weak attraction between the fibers, see SI). In 
this case, whose analysis is detailed in the SI, we do not observe significant deviations 
from the results presented in the main text. It is also worth noticing that this perfectly 
rigid environment closely resembles artificial gels made of solid nano-wires [34] which 
possess a much higher Young modulus and have been found to be optimal media for 
gel electrophoresis experiments. 

The external field is modelled as a force f acting on each bead forming the poly¬ 
mers. Assuming that in physiological conditions half of the charges from the phosphate 
groups are screened by counter-ions [38], one can think that each bead (cr = 2.5 
nm ~ 8 bp) contains a total charge of = 16qej2, where qe is the electron 
charge. Within this assumption we can map the external force applied onto each bead 
to an effective electric field E = f jqt- Although this mapping is a crude approx¬ 
imation of the Coulomb interaction between the charged DNA, the ions in solution 
and the applied electric field, we find that we can recover a “weak field” behaviour 
of the knotted samples, i.e. linear increase of the speed as a function of their ACN, 
up to ~ 50 Vjcm, which is roughly comparable with the field intensity used in 
experiments. In this work we used field intensities in the range from E = 1.25 
V jcm \.o E = 625 Vjcm. 

The average crossing numbers used in this work have been obtained from Ref. [39], 
where the authors computed the ACN corresponding to Mobius energy minimising 
knotted configurations. The thermally averaged ACN of the samples used in this work 
has been computed from equilibrated configurations and has been found to be in a 
one-to-one correspondence to the values in Ref. [39] (not shown), confirming the linear 
relationship between the ACN of ideal and thermally agitated configurations [16]. 

The hydrodynamics is here considered only implicitly, as it is customary for Brow¬ 
nian dynamics simulations. This means that the polymers do not feel one another 
via hydrodynamical interactions but are subject to thermal fluctuations due to a sur¬ 
rounding bath at fixed temperature T (see SI for more details). 

The simulation timescale is given in terms of the Brownian time, which corresponds 
to the time taken by a bead of size a to diffuse its own size, i.e. tbv = jE^, 
where = ksTj^ = ksT {S7Trjsol^)~^ is the diffusion coefficient of one 

bead and rjsol — 10 cP the solution (water) viscosity. From this we obtain 
TBr = STTrjsolcr^/^ bT ~ 40ns. 


ACKNOWLEDGMENTS. DMi acknowledges the support from the Complexity Sci¬ 
ence Doctoral Training Centre at the University of Warwick with funding provided 
by the EPSRC (EP/E501311). The computing facilities were provided by the Centre 
for Scientific Computing of the University of Warwick with support from the Sci¬ 
ence Research Investment Fund. DMa thanks EPSRC grant EP/1034661/1 for sup¬ 
port. EO acknowledges support from the Italian Ministry of Education grant PRIN 
No. 2010HXAW77. 

Appendix: Supplementary Information 


Model and computational Details. Nicked DNA knotted loops 
are modelled as coarse grained bead-and-spring circular chains 
with a knot manually tied in. Let Vi and di^j = Vj — Vi be 


respectively the position of the center of the z-th bead and the 
vector of length dij between beads i and j. 

The connectivity of the chain is treated within the finitely 
extensible non-linear elastic model [35] with potential energy, 


k 2 

UpENEii, z + 1) = ~ 



I di^i-\-\ 

V Ro 


2 


for di^ipi < Rq and UpENEidE + 1) = oo, otherwise; here we 
chose Rq — 1.6 (7 and k = 30 e/cr^ and the thermal energy 
kpT is set to e. The bending rigidity of the chain is captured 
with a standard Kratky-Porod potential, 


f/t,(z, z + 1, z + 2) 


kpTlx 

2(7 


di^ipi • <ii+i,z+2 

di^i-\-idipi^iP2 


where Ik = 40(7 ~ lOOnm is the known Kuhn length of uncon¬ 
strained DNA. The steric interaction between beads belonging 
to the polymers or the gel is taken into account by a truncated 
and shifted Lennard-Jones potential 


ULj{i,j) = 46 




6 

+ 1/4 


e(2^'^ac 


dij). 


where 0(x) is the Heaviside function and (7c can be, depend¬ 
ing on the case, (7c = cr (diameter of the chain beads) and 
o'c = (o'g E a)12 — 5.5 a where ag is the diameter of the beads 
forming the gel. Notice that steric interactions between beads 
of the gel are excluded from the computation since we consider 
the gel as a static mesh. This mimics artificial gels made of 
solid nano-wires, which have Young modulus 10^ times larger 
than that of agarose [34]. 

Denoting by U the total potential energy, the dynamic of the 
beads forming the rings is described by the following Langevin 
equation: 

mvi =— VU E T] [1] 

where ^ is the friction coefficient and r] is the stochastic delta- 
correlated noise. The variance of each Cartesian component 
of the noise, a^ satisfies the usual fluctuation dissipation rela¬ 
tionship ap = 2^kBT. 

As customary [35], we set m/^ = tlj, with tlj = a^/mje = 
ay/m/k^ and the Brownian time tbv — a j Dh with Db = 
kpT/^ is the characteristic simulation time step. From the 
Stokes friction coefficient of spherical beads of diameters a we 
have: ^ = 37r?7soZcr where 77 ^ 0 z is the solution viscosity. By using 
the nominal water viscosity, rjsoi = 1 cP and setting T = 300 
K and a = 2.5 nm^ one has rpj = rpr = SttPsoICt ^= 37 
ns. Note that, being the gel structure static for all times, no 
evolution equations for the beads belonging to the mesh are 
considered. The numerical integration of Eq. [ 1 ] is performed 
by using a standard velocity-Verlet algorithm with time step 
At = O.Olr^r ~ 0.4 ns and is implemented in the LAMMPS 
engine. 

The whole system is contained within a box of linear dimen¬ 
sion L = 320(7 with periodic boundary conditions in all three 
directions. The irregularities of a physical gel are created by 
starting from a regular cubic mesh and then by halving some 
of the edges randomly, with probability p = 0.4. The lattice 
spacing is set to I = 80(7 = 200 nm^ which is typical for an 
agarose gel at 5% [33], and the edges of the mesh are discre¬ 
tised with beads of size cr^ = lOcr ~ 25 nm, again compatible 
with the size of the agarose bundles [33, 37]. 

By considering the beads as cylinders with height and di¬ 
ameter equal to cr, the system volume fraction results = 
TvNMa^/4{L^ — Vgei) — 1.3 • 10“^, much smaller than the 
value at which the chains start to overlap {(fP ~ 2 • 10“^). In 
other words, the systems are in the dilute limit and interchain 
interactions occur rarely in the simulations. 
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Fig. SI. Panels (a) and (b): trajectory and squared radius of gyration vs time for some knotted DNA loops subject to weak (50 Vjcm) (a) and moderate (150 
Vjcm) (b) electric fields and moving through an irregular gel, i.e. with dangling ends. Panels (c) and (d) refer to the regular gel case with same lattice spacing. Note 
that for high fields and in gel with dangling ends the knotted loops are trapped (zero velocity) for long time in states with large radius of gyration, suggesting an entangled or 
“impalement” state, while in a regular gel the difference between the trajectories at weak and moderate fields is very little, (e) The velocity of knotted polymers in regular gels, 
i.e. without dangling ends, displays a monotonic behaviour as a function of the average crossing number. For dense enough gels {I = 200 nm), the contribution coming 
from the interaction with the structure is observed to be sufficient to recover the linear electrophoretic separation as a function of the average crossing number at weak fields 
{E = 1.25 Vjcm). This is instead not observed in the case of sparser gels {I = 500 nm), where we recover no electrophoretic separation (f). We observed a deviation 
from the linear trend only at stronger fields and denser gels {E > 12.5 Vjcm, I = 200 nm), although we never observed a field-inversion or non-monotonic behaviour 
(g) 


The starting point of the various simulations is a configura¬ 
tion in which simple knotted configurations are placed outside 
the box. We then slowly pull the rings inside the gel structure, 
avoiding any impalement. Finally, we place the boundaries of 
the simulation box so as to match the boundaries of the gel 
structure and impose periodicity along all the three coordinate 
directions. To allow the equilibration of the system from this 
initial state, we first run 1 • 10^ TBr time steps, which are dis¬ 
regarded before the external force is switched on. The mean 
and standard error of various properties, as the mean displace¬ 
ment along z, are calculated by averaging over the rings in the 
system and over different runs starting from different (equili¬ 
brated) initial conditions. 

The magnitude of the force / acting on the beads is ex¬ 
pressed as a multiple of the system units e/cr ~ 1.6pN. The 
force acting on the beads can therefore be expressed as the 
resultant of an electric field E = 1.6/qb pN ~ 12.5 kVjcm. 
In this work we have used forces / in the range between 


e/a and 5 e/a, which can then be mapped to electric 

fields roughly between 1.25 V/cm and 625 V/cm, compati¬ 
ble with the values used in standard DNA gel electrophore¬ 
sis [40, 41, 12, 4]. 

We compute the mobility of knotted polymers by measur¬ 
ing the centre of mass displacement along the direction of the 
applied field. In Fig. SI we show single example trajectories 
of the knots and their relative change of radius of gyration 
during in silico gel electrophoresis experiments. In Fig. SI, 
we also compare trajectories and speed of knots travelling in 
regular and irregular gels. As one can notice, while in irregu¬ 
lar gels (Fig. Sl(a)-(b)) the motion of the knots under weak 
(50 V/cm) and moderate (150 V/cm) fields displays a rather 
different behaviour, in regular gels (Fig. Sl(c)-(d)) knots that 
are driven by the same external fields display similar trajecto¬ 
ries. In particular, we do not observe a field inversion at any 
field or gel density (Fig. Sl(e)-(g)). 
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Fig. S2. (a)-(g) Waiting time distribution -P(Tdis) for different knot types and in units of simulation Brownian time Tsr- can notice that the distribution shows 

to peaks at short and large times which correspond to “transverse” or “head on” interactions, respectively, (h) Total amount of time spent in an entangled state as a function 
of the ACN. This quantity reflects the fact that slower knots spend more time in entangled states (compare with Fig.2 in the main text). 


The role of Hydrodynamics. The Brownian dynamics scheme 
allows us to introduce an implicit solvent, i.e. there is no 
explicit hydrodynamics taken into account during the simula¬ 
tions. We made this choice assuming that in the case of dense 
gels, the contribution coming from hydrodynamics would be 
negligible compared to the interactions of the polymers with 
the environment. This is conhrmed in Fig. SI where we show 
the linearly increasing speed of knots travelling through a reg¬ 
ular gel at weak and strong fields. In both cases, even with¬ 
out the presence of hydrodynamics or rescaling, we recover 
the expected linear relationship between speed and ACN. In 
Ref. [15] the authors rescaled the hnal velocities of the knots 
by using the Kirkwood-Risenman formula, thereby accounting 
for hydrodynamics. In this work we observed that such rescal¬ 


ing does not affect the hnal functional form of the results, con¬ 
sequently we did not introduce such rescaling. We argue that 
the contribution from hydrodynamics becomes relevant when 
sparser gels are considered. In this case we observed (data not 
shown) that the rescaling operated in Ref. [15] is necessary 
to recover the experimental observations as no electrophoretic 
separation is observed (see Fig. Sl(f)). 

Waiting Times. In order to quantify the time taken by the 
rings to disentangle from the impalements, we analyse the sin¬ 
gle trajectories (reported for example in Fig. SI) and compute 
the amount of time each knot is moving slower than it would 
if it were free, i.e. Vfree — f/C- From this we can quantify 
the total amount of time each knot is stuck in an entangled 
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x = r/l 

Fig. S3. Probability for a circle of radius r drawn randomly on a square lattice with lattice spacing I to encircle a lattice vertex {Py - black circles) and to overlap over a 
lattice edge (Pe “ gi'ay squares). The two quantities reach 1 for r = Z/2 (Pe) and for r = Ij\/2 ~ 0.7 {Py). The functions and parameters used to fit the data are also 
shown. 


state. This is found to be a non-monotonous function of the 
ACN (see Fig. S2(h)) and is in line with the non-monotonic 
average speed of the knots and the relative separation in gel 
electrophoresis experiments as a function of their ACN. 

Each time a knot is stalled, it takes a certain amount of 
time to re-establish its motion. We report the distribution of 
these “waiting times” (rdis) in Fig. S2(a)-(g) for several knots. 
As one can notice, the distribution is bimodal, i.e. shows two 
distinct peaks at short and long times. We interpret the peak 
at short time-scales as transversal impalements or crawling 
around the gel structure, while the peak at long times as the 
fact that sometimes knots interact via an “head-on” collision 
with the gel open strands. This immobilises the knots for long 
time, as the only way to re-establish the motion is by moving 
against the field for a length of at least half lattice spacing. 

Biased Continuous-Time Random Walk Model. Here we de¬ 
tail the biased continuous time random walk (CTRW) model, 
whose results are reported in the main text. 

We start by noticing that the simulations suggest the exis¬ 
tence of two important time-scales in the dynamics: r/ and 
Tdis, respectively regulating the frequency of entangling events 
and the time required to disentangle from the impalements 
with the dangling ends. Firstly, we address the situation in 
which Tdis ^ 1. In this case the knots disentangle quickly, and 
we can think of the whole process as random walks biased in 
the opposite (the DNA is negatively charge) direction of the 
field. Due to topological interactions with the gel, the biased 
random walkers eventually stop with a rate Ae — We 

expect that this rate should be described by the probability 
Pe{r/l) of a circle of radius r to overlap onto an edge of 2D 
square mesh with lattice spacing /, i.e. 

where x(JC) = Rg(JC)/l. Note that this rate is normalised 
with respect to the empirical value observed in the Brownian 
dynamics simulations for the unknot (Oi) at weak fields. 

We can obtain the value of Pe{x) as a function of the ratio 
Rg{lC)/l^ with I the gel lattice spacing from Fig. S3, where we 
simulated the deposition of a circle of radius r onto a 2D mesh 
with lattice spacing I and calculated the probability of over¬ 
lapping over a strand {Pe{x)) or to encircle a vertex (P^;(x)) 
of the grid. 

The knots moving in the gel are therefore viewed as random 
walkers on a lattice. They can jump one site to their right 
(broadly this corresponds to one unit cell of the gel), or are 


required to stop for one time-step (because Tdis ^ 1) with rate 
Ae- This is equivalent to a Poisson process whose probability 
of getting entangled after t time-steps is given by: 

Pent{x{ACN),t) = 1 - [1 - Xe{Rg{IC))Y ~ 1 - exp [-Xe{IC)t], 

[3] 

where we stress the dependence of Ae on the knot type via its 
size Rg. The average time between interactions is given by 
(t/) A“^ and the average number of stops on a time t is 

{ue) = Aet. The relative separation between two knots can be 
expressed as 

Ax(/Ci, /C2) = Vf ree (tf ree ree (^2)) 

— Vfreet(Xe(JC2) “ Ae(/Cl)) [4] 

where Vfree is the velocity of the knots when no interactions 
occur (when hydrodynamics is negligible this is the same for 
all knots), t is the observation time and t/ree — t — Xet. From 
this is clear that the separation between knots JC 2 and JCi is 
proportional to the difference of their entanglement rates, as 
expected. Given that our model relies on the hitting prob¬ 
ability between a knot and its environment, it is clear that 
if the knot size is either much larger or much smaller than 
the gel lattice spacing, the electrophoretic separation is not 
as good as in the case the two are comparable. Within our 
simple model we recover good linear relationship within the 
range Rg/l E [10“^, 10] meaning that short samples (or large 
gel pores compared to the gyration radius of the knots) work 
better for the electrophoretic separation at weak fields rather 
than small pores (or very large Rg) [29] (see Fig. S4(a)). This 
is compatible with the fact that is often impractical to per¬ 
form gel electrophoresis with very large samples. In fact, for 
other cases outside this range, we observe essentially no spatial 
separation of the knots (see Fig. S4(a)). 

We now introduce a topology dependent disentanglement 
time: every time that the Random Walker undergoes a “stop”, 
we require that it has to wait an amount of time that increases 
proportionally to its average entanglement number (see (tt) in 
Fig. 3 of main text). We pick the disentanglement time from 
an exponential distribution with characteristic time given by 
the empirical one observed in the MD simulations (see Fig. 4 
in the main text), i.e. 

tw — In t/ Adis [5] 

where r is a random number and A^-J = Tdis- As observed in 
the main text, the disentanglement time Tdis is expected to be 
a function of (tt), i.e. to increase linearly with the ACN. We 
therefore use: 

rdis(ACN) = A{7t) + Tdis(Oi) [6] 
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Fig. S4. Average relative separation over time (in units of lattice spacing over time) of the knots as a function of the ACN obtained from the simple CTRW model described 
in the text, (a) For weak fields (for which the disentanglement time is quick T^is 1) the linear relationship is re-established for Rg ~ I, while in the two limiting cases 

Rg ^ I and Rg Z we do not observe a clear separation, as expected, (b) The non-monotonic behaviour of the separation is re-established when a topology-dependent 
disentanglement time is introduced in the model. The typical electrophoretic arc crucially depends on the ratio Rg(ACN) /I. (See text for details). 






ACN 



Fig. S5. (a) Picture of the gel with flexible dangling ends (light grey) anchored to the rigid structure (dark grey). (b)-(c) Average velocity of some knots under weak 
(50 Vjcm) and moderate (150 Vjcm) field, respectively, (d) Snapshots (not to scale) of entanglement events of 3i (top) and 9i (bottom) knots corresponding to 
stalling points in their trajectories and large values of Rg. 


where A = 5 rsr is an empirical parameter fitted from the 
simulations and Tdzs(Oi) = SOr^r is the empirical average dis¬ 
entanglement time of the unknot. In addition, long waiting 
times which are typical of head-on impalements are taken into 
account by adding a small probability {q ^ 0.05, compatible 
with the bimodal distribution in Fig. 4 of the main text) and 
conditioned to the fact that an entanglement event happened. 
For this we set Tiong = A(7r) + 250. So the average disentangle¬ 
ment time from a head-on impalement is still knot-dependent 
but is much longer than any characteristic disentanglement 
time for other types of entanglements. 

Within this model we recover the electrophoretic arc in 
Fig. S4(b). In addition, the shape of the arc can be hnely 
tuned to obtain different depths of the arc and different slow¬ 
est knots. This is mainly controlled via the ratio Rg/l (we 
here assume that the external held is neither too weak or too 
strong). In Fig. 5 of the main text we show that within this 
model we can accurately recover the electrophoretic speed of 
the knots observed in the MD simulations by shifting the curve 
obtained with the CTRW model setting R^(Oi)// = 0.5 (as in 
the MD simulation) to the observed speed of the unknot and 
by rescaling it by the free speed vfree at which the knots travel 


between two planes of the gel, i.e. two lattice sites in the 
CTRW model. It is worth noting also that this model allows 
us to predict the mobility of knots which were not probed by 
Brownian dynamics, by making use of the empirical functional 
forms of (tt) and Rg derived from these simulations. It is also 
worth stressing that the only real free parameter of this model 
is the lattice spacing /, while the value of R^(Oi), Tdis(Oi) and 
Ae(Oi) are input parameters needed to calibrate the model to 
a target case. 

Gel with flexible Dangling Ends. While the assumption that 
the gel is perfectly rigid is only an approximation for agarose 
gels, several experiments reported values for the persistence 
length of single agarose hbers in the range of 2 — 10 nm; in 
particular Ref [37] reports /agarose = 9 nm (the same work 
also presents evidence of the presence of dangling ends in 
the gel). Furthermore, each agarose bundle (or “fibril”) is 


^The simulations in Fig. S5 consider flexible bundles with fixed angles at crosslink, i.e. the clamping 
at the crosslinks is infinitely strong. For the case of the unknot, we checked that replacing this right 
angle constraint with a soft angular potential keeps the dynamics qualitatively unchanged (we still 
observe entanglements or impalements). 
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formed by a number n/ of fibers, with n/ ranging between 10 
and 20 [33, 37], and has a diameter cr/ ~ 10 — 20 nm [42]. 
The persistence length of a full fibril might depend linearly or 
quadratically on n/ (depending on the structure of the bundle 
and strength of interfiber interactions, see Ref. [43]). There¬ 
fore, a conservative estimate of the persistence length of an 
agarose bundle is attained if a linear scaling is assumed (cor¬ 
responding to weak attractive interactions between the fibers, 
see Ref. [43]), leading to Agarose Idna — 50 nm. 

In order to probe the effect of fibers’ flexibility on the en¬ 
tanglement/disentanglement process experienced by the knots 
travelling through the gel, we therefore modify the gel struc¬ 
ture to accommodate flexible dangling ends (see Fig. S5(a)) 
described as semi-flexible polymers with persistence length 
equal to that of dsDNA, and anchored to cross-links that 
are still modelled as rigid and static. This choice corre¬ 
sponds to flexible agarose bundles clamped at right angles at 
crosslinks and avoids the un-physical scenario in which the 
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